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SIMIAET 


A formula, Based on a semlrational analysis, is presented for 
estimating the torsional stiffness necessary to prevent flutter of a 
swepthack or unswept uniform wing that attains supersonic speeds . 

■ Eestilts of missile flints at speeds up to Mach number 1.4- demonstrate 

the usefulness of the formula. 


INTROIXJCTION 


FallTAres probably due to flutter were encountered in NACA flif:Ji'b ‘ 
tests of several rocket-powered, drag-research, missiles that were 
intended to attain Meoh numbers of about 3..k. The wing failures of 
these missiles led to the development of a simple, semiratlonal 
torsional stiffness criterion for preventing flutter of uniform, 
sweptbaok or unswept missile wings that attain supersonic speeds. 
Missiles that failed were redesigned in accordance with this stiffness 
criterion and proved to be safe in flight. 


TOESIONAL STIFFNESS CSITEBION 


On the basis of the semiratlonal an^y sis presented in appendix A, 
the following formula is proposed for estimating the torsional stiffness 
necessary to prevent flutter of a unifown sweptbaok or unswept wing 
that attains supersonic speeds; 
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where (fig. 1) 


GU torsional stiffness (ratio of torque to twist per unit length) 
of section normal to leading edge^ poiind -inches^ 

L length, inches 

c chord, normal to leading edfe, inches 

4 Distance of center of gravUy "behind quarter-chord position 

Chord 

Equation (l) may be considered as probably most reliable for wings 
having the following characteristics: 


(a) Low ratio of bending frequency to torsional frequency: 

1 

(b) High relative density: > 10 (see appendix B) 


(c) Center of gravity ahead of midchord position. However, for 
wings that do not quite satisfy those co:idltions, the criterion may 
be used as a design guide until more experimental find theoretical 
information becones available. 


The derivation of equation (l) was mdo for standard sea-level 
atmospheric* conditions , application of the formula to high-altitude 
conditions is probably conservative. 


Divergence of Unewept vrinGs 

Failure by divergence ratlier than by flutter may occur in wings 
without sweepback. Let e be defined by; 

Dieteinoe of shear center behind qxaar ter -chord position 

Chord 


then the larger of d and e should be used in oqmtion (l) in 
order to guard against the possibility of divergenco as well as flutter 
of unswept wings. 
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FLIGHT lESTS 


A ntmiber of roclret -powered iniasilas with xmlfoiTa wooden vringa, 
sweptback' and laiswept, have been flown by the KA.CA in the covirse of an 
investigation of drag at speeds up to ^^ach niaaber lA. Gome of the 
missiles lost theix’ wings in fli^t; subsequent models of these missiles 
were flown successfully'’ af-fcer the wings •were reiatforced wi’tli almninum 
sheet bonded to -fche upper and lower surfaces in order to increase ■fche 
torsional stiffness sufficiently to satisfy the criterion presented in 
the present paper. 

A history of the fli^t experience with missiles is summarized 
in figure 2, which compares ■the actizal wing stiffnesses (measured or 
calcula’ted) with the stiffnesses required tty prevent flutter according 
to equation (l) . The data for figure 2 are sho^m in table I. It is 
to be no’tod "^at all missile \rinfjs 2 wi-fch •torsional stiffnesses •tliat fall 
above the straight-line plot of -the stiffneos cri’torion did not fail 
in flight. The presence below the line of the two points representing 
missiles that did not fall indicates some conser'vatism of the 
criterion. 


CONCniDII'fG KEMAEK 


The usefulness of the torsional stiffness criterion presen-ted 
has been demonstra’ted by the rooilLta of a limited number of flight 
■tests of missiles with uniform^ sweptback and unswept wings. Ho’ifover, 
the criterion should be regarded as subject to modification or replacement 
as experlmen'tal and -theoi’otical da-ta in greater quantity and at speeds 
higher than ^■Iach number l.it bocomo a-vailable. 


Langley Memorial Aeronautical Twaboi’atory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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APPEMDIX A 

DERIVATION OP TORSIONAL STIFFNESS CRITERION 


The following analysis refers throu£^iout to unswept wings. 
However, oonsidora'ble unpublished NACA data, as well as the data 
of reference 1, indicate that a wing of given L and c (fig. 1(a)) 
has a higher flutter speed if it is sweptback than if it is unswept. 
Hence, the stiffness criterion developed should be conservative when 
applied to sweptbacl: wings . 

Flutter at low subsonic speeds . - In reference 2, The odor sen and 
Garrick present the following empirical formula for the flutter at 
low subsonic speed of a two -dimjsnsicnal wing (fig. 1(c)): 


^ 

b%, y ** ^ + a + 


The formula is stated to be reasonably good for wings having small 

and small k . The following theoretical formxila for the divergence 
speed of a two-dimensional wing is also given in reference 2J 



(A2) 


Using strip theory, several authors (references 3, and. 5) have 
given as the divergence speed of an unswept uniform three-dimensional 
wing (fig. 1) 



(A3) 


Equations (Al) and (A2) differ only in that the term '^ + a in the 

divergence eqtiation is replaced by ± + a + 2^ in the flutter eqmtion. 

Then, by analogy with eqiiations (Al) and (A2) , equation (A3) can be 
modified to give the flutter speed of a iiniform three-dimensional wing 
by replacing o by d; thus. 
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The value of to he used in e(juatlon (A4) is the two -dinonsional 

value 2rt multiplied hy an aspect-ratio cori’ection. In caloxilating 
divergence speeds hy strip theory, Shomiok (reference 6) mal’.es the 
approximation 


or 



(A5) 


(For the case of a uniform wing with Ir/c = 3*1^, this assun^ition 
gives for the divergence speed computed hy strip theory, eq^uation (A3) , 
a value that differs hy less tlian 1 percent from the exact result 
calculated hy lifting-line theory hy Hildehrand and Reiesner (reference 5) • ) 
Use of equation (A?) in the flutter oq,uatlon (A4) gives for the req.uired 
torsional stiffness to prevent flutter at low suhsonic speeds 


GJ « 



(A6) 


Flutter at hifh suhsonic sneeds .- Eq.uation (A<5) may he extended 
to speeds up to about M = 0.75 hy using the Glauert-Prandtl 


con^iressihillty correction on 



•thus, in eq.Tiation (A6) let 
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This use of ■Hi© Glauert-Prandtl correction is sucgested "by Garrlci: 
(reference 7) for vings with small suad small k . From 

Garrick's stvuiy of numerical flutter calculations in reference 7j it 
may 1)© concliaded that for other types of iiJlngs this corroction is 
conservative . 

Flutter at transonic and suTaersonic sueeds . - From the studies 
of Garrick and Fuhinow (reference 8) on supersonic flutter, the 
following conclusion may be droim; For wings having low 

and low it , and having the center of gravity eihoad of the midchord 
position, ihe transonic range appears critical for bending-torsion 
flutter. That is, if the wing passes thi'ough the transonic range 
without fluttering, it will probably not flutter at hl^er speeds. 
Furthermore, even if the conditions specified on ccji/ttia,, and 

the position of the center of gravity are not wholly fulfilled, it is 
probable that if the wing passes safe3.y throu;^ transonic speeds, it ‘ 
will not flutter until a Mach number considerably hlghor thaii 1 is 
attained . 

With these considerations in mind, it appears that a procedure to 
prevent flutter of a large class of suTJersonic missiles is to design 
against ■transonic flut'fcor. Since -there is very little transonic 
aerodynamic information available, the a^-bhod 'to be used is to extend 
the form of eq.ua tion (a 6) . For the purpose of the present analysis 
it will be assumed that equation (A7) holds up to M = 0.75, and 
•that between M = 0.75 and M = •! (fig. 3), 


. ..-,^1 

>^/oo /1 - (0.75)"^ 


= 2t:(1.51) 


The design value of v^ wi?.l be taljen as the velocity of sound. Then, 
substituting in the design formula (a 6) the values 


Vf = 1120 ft /sec 


(Velocity of sound at 
sea level) 


p s 0,002378 lb-sec2/ft^ 


(Density at sea level) 
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gives 



« 2jt(l-5l) 



2 

■where L and o are in inches, and GJ in pound **inohes • 
Bounding off the value of the Qons’fcant gives as ■tiis final design 
formula 


GJ = ItO 


L3c2d> 



HICVJ HICVJ 
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APPENDIX B 
SYMBOLS 


c 

d 


e 


L 


Vf 

GJ 


M 

P 

Bcl 


A 


chord normal to leading odee (See fig. 1.) 

Distance of oentor of gravity 'behind oTiarter-chord -position 

Chord 

PistEince of shear center "behind quarter-chord position 

Chord 

length along leading edcp (Gee fig. 1.) 
divergence speed 
flutter speed 

torsional stiffness, ratio of torque to twi.st per unit 
length 

Mach number 

air density 

lift -curve slope, f ini to wing 
lift -curve slope, infinite wing 

angle of sweepbacl: 


The followisag symbols and their definitions eire essentially those 
of Theodorsen and Garriclv, reference 2: 

b half chord, used as reference vuilt length 

Distance of shear center behind qiteirter -chord pdsitlon 

Half -chord 

(See fig. 1.); ^+a=2o 

+ a + x_ P3. stance of center of Gravity behind, quarter -chord position 
^ Half -chord 

(See fig. 1.); ■^ + a + = 2d 



Maas radius of gyration referred to ahe&r center 

Half -chord, 

ratio of mss of cylinder of air of diameter eq.ual to 
chord of wine to mss of wing, both talien for eq.un.’. 
length along span 

angular frequency of uncoupled toreional vibration about 
shear center 

angular frequency of uncoupled bending vibration 
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lABLE I 


DATA PLO!HED IN FIGURE 2 


Missile 

A 

(deg) 

L 

(in.) 

c 

(in.) 

d 

GJ 

(kip -in. 2) 

Fliekt 

reeiilt 

(a) 

1 

0 

10-37 

9.63 

0.18 

67.3 

N 

2 

0 

12.92 

8.12 

.18 

52.2 

N 

3 

0 

14-03 

6.88 

.18 

21.9 

F 

k 

0 

14.63 

6.88 

.18 

145.0 

N 

5 

0 

10.37 

9.63 

.25 

7li6 

N 

6 

34 

8.63 

10.67 

.18 

106.0 

N 

7 

34 

12.51 

7.90 

.18 

31.8 

N 

8 

34 

15.25 

6.56 

.18 

24.6 

N 

9 

34 

15.25 

6.56 

.18 

109.7 

N 

10 

34 

17.63 

5.69 

.18 

11.4 

F 

11 

34 

17.63 

5.69 

.18 

95.2 

N 

12 

45 

10.11 

9.10 

.18 i 

55.8 

N 

13 

45 

14.75 

6.88 

.18 

18.3 

N ■ 

14 

45 

18.28 

5.74 

.18 i 

13.0 

F 

15 

45 

20.63 

4.88 

.18 

6.8 

F 

16 

45 i 

14.75 1 

6.88 i 

.25 ! 

20.4 

F 

17 

45 

14.75 

6.88 

.25 

147.6 

N 

18 

52 

11.61 j 

7.94 

.18 

32.5 

N 

19 

52 

16.88 

5.88 

.18 

9.7 

F 

20 

52 

16.88 

5.88 

.18 

91.6 

N 

21 

52 

21.00 

5.00 

.18 

7.6 

E 


®N - No fallxjre 
F - Failure 


MTIOKAL ADVISORY 
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( 0 ) Uniform sweptbock wing. 



(b) Notation of this report* 



Figure 1.- Symbols for wing dimensions. 
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Figure 2.- Comparison of torsional stiffness criterion 
with flight test experience. 




Lift curve slope, 


14 


NACA RM No. L7G02 



Mach number 


Figure 3.- Assumed variation of lift curve slope 
with Mach number for purposes of 
flutter analysis. 
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